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Abstract 

Aims. We study the differences between non-magnetic and magnetic regions in the flow and thermal structure of the upper solar 
photosphere. 

Methods. Radiative MHD simulations representing a quiet region and a plage region, respectively, which extend into the layers around 
the temperature minimum, are analyzed. 

Results. The flow structure in the upper photospheric layers of the two simulations is considerably different: the non-magnetic simu- 
lation is dominated by a pattern of moving shock fronts while the magnetic simulation shows vertically extended vortices associated 
with magnetic flux concentrations. Both kinds of structures induce substantial local heating. The resulting average temperature pro- 
files are characterized by a steep rise above the temperature minimum due to shock heating in the non-magnetic case and by a flat 
photospheric temperature gradient mainly caused by Ohmic dissipation in the magnetic run. 

Conclusions. Shocks in the quiet Sun and vortices in the strongly magnetized regions represent the dominant flow structures in the 
layers around the temperature minimum. They are closely connected with dissipation processes providing localized heating. 
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1. Introduction 

Semi-empirical models of strongly magnetized regions in the so- 
lar atmosphere (such as network and plage areas) indicate in- 
creased t emperatures relative to the 'quiet' Sun (e.g. |Fontenla| 
et al.|2006| ). This is consistent with the fact that individual mag- 



netic structures are hotter than their surroundings in the upper 
layers of the photosphe re and in the chromosphere (e.g.|Solanki 
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|1993[|Lagg et al.|2010| ), but the nature of the concomitant heat- 
ing process is still under debate. One possibility is the dissipation 
of mechanical energy in the form of flows or waves, another is 
resistive dissipation of magnetic energy, for instance in current 
sheets. 

The dynamics in these atmospheric layers is complex and di- 
verse. For the weakly magnetized internetwork regions, numer- 
ical simulations and observations (e.g. Carlsson & Stein 1992 



[19971 |Wedemeyer et al.|2004[ |Wedemeyer-Bohm|2010f suggest 
that shocks driven by overshooting convection play an impor- 
tant part above the temperature-minimum region. In addition, 
swirling m otions have been detected in observations of the chro - 
mosphere ( Wedemeyer-Bo hm & Rouppe van der Voort|2009|), 
the transition region ( |Curdt et al.|2012| ), and the corona ( Zhang 
|& Liu|2011| ). 

Small-scale vortices occuring in simulations of non- 
magnetic or weakly magnetized regions of the lower atmosphere 



and uppermost convection zone were studied by Moll et al. 
( 2011 ) and Kitiash vili et al.| ( [20TT| ). Simulations comprising suf- 
ficient magnetic flux to represent network patches or plage re- 
gions exhibit vertically orientated vortices associated with mag- 
netic flux concentrations in the upper photosphere and lower 
chromosphere (|V5 gler 2004 ; Shelya get al.|201 l[|Carlsson et al 



EOlOUKitiashvili et al. 



20T2l|Steiner & Rezaei|20T2] )7 



In this paper, we analyze simulations extending higher into 
the lower chromosphere for a detailed study of the effect of a 



magnetic field on the flow structure as well as on the heating 
processes in the atmosphere up to a few hundred km above the 
average height of the temperature minimum. We compare a sim- 
ulation without magnetic field as an extreme case of 'quiet Sun' 
with a simulation having a mean vertical field of 200 G, which 
is considered to represent the conditions in a plage region. 



2. Methods 

2.1. Simulations 

We analyzed simulations carried out with the MURaM code (for a 
detailed description, see Vogler 2003 ; Vogl er et al.|20 05 ), which 
treats the equations of compressible MHD together with the ra- 
diative transfer equation solved along 12 rays per grid cell. The 
equation of state incorporates the effect of partial ionization for 
1 1 elements. The numerical scheme uses 4th-order centered spa- 
tial differences and a 4th-order Runge-Kutta method for the time 
stepping. Both the top and bottom boundaries permit free in- and 
outflow of fluid, the side boundaries are periodic. The magnetic 
field at the top and bottom boundaries is assumed to be vertical. 

Here we consider two simulation runs, one non-magnetic 
and one with an average unipolar vertical magnetic field of 
200 G. Both runs used a computational box with a horizon- 
tal area of 6 x 6 Mm 2 and a height of 1.68 Mm, ranging from 
about 900 km below the average height of the optical sur- 
face to roughly 800 km above. The size of the grid cells was 
20.8 x 20.8 x 14.0 km 3 . The simulations were run for several 
hours solar time to reach a statistically stationary state. For the 
statistical analyses, we used 4 snapshots from each run. 
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Figure 1. Volume rendering of the swirling strength, A^, in the non-magnetic case (left) and the magnetic case (right) from two 
viewing angles. Shown is the upper half (6 x 6 x 0.8 Mm 3 ) of the computational domain, the bottom plane corresponding roughly 
to the optical surface. The color scale covers the range < A c { < 0.027 s" 1 ; higher values of A c { are saturated. 



2.2. Vortex identification 

The vorticity is not the optimal quantity to identify swirling 
flows. For instance, a purely bidirectional shear flow without any 
rotational component nevertheless has a non- vanishing vorticity. 
For the detection of vortices, i.e., fluid elements rotating about 
a local, p ossibly mov i ng axi s, we followed the procedure em- 
ployed in Moll et al. ( 2011| ): vortices (swirling flows) are de- 
fined as regions where the velocity gradient tensor has a pair of 
complex conjugate eigenvalues ( [Zhou et al.|1999| ). The strength 
of the vortical motion is determined by the swirling strength, 
A C {, which is the magnitude of the unsigned imaginary part of 
the complex eigenvalues. In the case of rigid rotation, we have 
A C { = 27r/r s , where r s is the rotation (or swirling) period. The 
direction of a vortex is defined to be along the eigenvector cor- 
responding to the real eigenvalue. The inclination angle of this 
vector with respect to the vertical direction is denoted with l. 



3. Results 

3.1. Vortex properties 

Figure [T] shows regions of high swirling strength in snapshots 
from the non-magnetic and magnetic simulations, respectively. 
The most striking difference between the two cases is the domi- 
nance of tall, vertically orientated vortices in the magnetic case, 
which extend over the entire height of the simulated photo- 
sphere. In the non-magnetic case, the vortex features do not pro- 
trude far above the optical surface, bending over to form mainly 
horizontally orientated loops. Vortex features near the top of 
the computational domain are associated with shock fronts (see 
Sec. |3.2| ). The difference between the two simulations is fur- 
ther illustrated in Fig. [2j which shows the location of vertical 
(i < 20°, left column) and horizontal (i > 70°, right column) 
vortices on horizontal planes. In the upper two rows, the vortex 
locations are compared to the distribution of vertical velocity at 
the average height of the optical surface (z = 0). Near the optical 
surface, both simulations show horizontal vortices at the edges 
of granules and vertical vortices in the intergranular downflows 
lanes. The comparison with the magnetic-field distribution in the 
third row indicates that the horizontal vortices near z = are a 
hydrodynamical phenomenon of the overturning motions at the 



borders of granules, which is essentially u nrelated to the inter- 
granular magnetic flux concentrations (cf. Steiner et al. 2010; 
Moll et al. 2011). In contrast to this, the vertical vortices are 
almost exclusively found in magnetic flux concentrations. This 
one-to-one relationship extends throughout the atmosphere, as 
shown in the bottom row of Fig.|2j which corresponds to a height 
of z = 644 km above the average optical surface. At this height, 
we find a multitude of vertical vortices in the (expanded) flux 
concentrations, while strong horizontal vortices are absent. 

The visual impression provided by Figs.[T]and 2]is confirmed 
quantitatively in Fig.|3j which shows height profiles of horizon- 
tally averaged vorticity, swirling strength, and the shear part of 
vorticity, oj^h = oj - 2A C {, where oj is the modulus of the vortic- 
ity vectorQln the magnetic simulation, all three quantities are 
suppressed in the deeper layers up to z - 200 km, owing to the 
effect of the Lorentz force on the convective motions. In the up- 
per layers, the situation is reversed and the swirling and shearing 
motions are enhanced in comparison to the non-magnetic case. 
We surmise that this is caused by the dynamical coupling of the 
vortex motion in the convective near- surface layers to the upper 
layers, which is mediated by the mainly vertical magnetic field. 

The drastic change in the properties of the vortex flows above 
about 150 km height is analyzed in more detail in Fig. [4] which 
shows height profiles of the fractions of horizontal area cov- 
ered by swirls of different strength and orientation. In the non- 
magnetic case (upper panel), horizontal vortices dominate at all 
heights, but the area fraction of strong vortices (swirling period 
r s < 2 min) drops drastically above the optical surface. For the 
vertical vortices, the decrease starts even below the optical sur- 
face. The increase of horizontal vortices above about 500 km 
is probably related to the pattern of shock waves developing in 
these layers (see Sect. |3.2| ). In the magnetic case (lower panel of 
Fig. [4]), the horizontal vortices behave very similar to those in 
the non-magnetic simulation, their area coverage dropping even 
faster with height. However, the area fraction of vertical vortices 
shows almost no drop above z - 0, followed by a strong increase 



1 For a purely rotational flow, we have oj = 2/L ci = oj r at the cen- 
ter of the vortex. Therefore, one can formally split the vorticity into a 
rotational part, oj r , and the remainder, oj sh = oj - oj r , which gives the 
contribution from locally non-rotating shear flow. 
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Figure 2. Distribution of vertical vortices (inclination i < 20°, left column) and horizontal vortices (i > 70°, right column) on 
horizontal planes. Pixels corresponding to swirling periods r s < 240 s are indicated in green color on cuts of the vertical velocity at 
z = (average height of the optical surface) for the non-magnetic case (top row) and for the magnetic case (second row). Similarly, 
the strong swirls are indicated on maps of the vertical magnetic field component at z = (third row) and at z = 644 km (bottom 
row) for the simulation with magnetic field. 

above 200 km, from which height on they dominate. Virtually Fig. [5] shows histograms of vortex properties at three height 
all swirls at large heights have vertical orientation (^<20°) in ranges, one in the deep layers of the computational domain 
the magnetic case. Even very strong vertical vortices with pe- (brown curves), one straddling the optical surface (green curves), 
riods below one minute reach appreciable area coverage above and one high in the photosphere (blue curves). While the vor- 
z - 400 km. tices in the non-magnetic case (dashed curves) and the magnetic 
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Figure 3. Horizontal averages of vorticity (top), swirling 
strength (middle), and shear part of the vorticity (bottom), as 
functions of height. Solid lines correspond to the magnetic case, 
dashed lines to the non-magnetic case. 



case (solid curves) do not strongly differ in the deep and the 
near-surface layers, they show strong deviations in the higher 
layer: in the magnetic run, the typical swirling periods are much 
smaller (i.e., the vortices are more vigorous) and the gas pressure 
is strongly reduced in comparison to the average gas pressure 
at the same height. The low pressure results from the location 
of the vortices in magnetic flux concentrations (cf. Fig. [2j bot- 
tom left panel), together with the effect of the centrifugal force. 
The temperatures of the weak swirls in the upper photosphere 
of the non-magnetic simulation are significantly higher than the 
average temperature, presumably due to their association with 
shock fronts. The vortices in the upper layers of the magnetic 
simulation are also hotter by 10-20% compared to their aver- 
age surroundings. The vertical velocities in vortices do not dif- 
fer strongly between the two simulations: in the deep and near- 
surface layers, vortices are predominantly located in downflows 
in both cases, while in the upper layer there is a slight preference 




Figure 4. Height profiles of the horizontal area fraction of grid 
cells with swirling period, r s = 2n/A C i, below given thresholds 
in the non-magnetic case (upper panel) and in the magnetic case 
(lower panel). Full lines correspond to vortices of all inclina- 
tions, while dotted and dashed lines separately represent hori- 
zontal and vertical swirls, respectively. 



for upflows in the non-magnetic case and an almost symmetric 
distribution around zero in the magnetic run. 

3.2. Dynamics of the upper photosphere 

High in the photosphere and above, the thermal and flow struc- 
ture is very different in the two cases. Figs. [6] and [7] show maps 
of various quantities on horizontal planes. In the non-magnetic 
case (Fig. [6]), there is a pattern of shock fronts with regions of 
strongly localized heatin g, leading to a thermal bifurcation into 
hot and cool regions (cf. Wedeme yer et al.|2004| ). The generally 
inclined and propagating shock fronts form a rapidly evolving 
filamentary pattern on horizontal cuts of temperature and diver- 
gence of the horizontal velocity. Some of the shocks are asso- 
ciated with vortices, which are preferentially horizontally orien- 
tated (cf. Fig. |4} upper panel). An example of a shock structure 
is given in Fig7[8j showing that the post-shock peak temperatures 
can reach values around 7000 K, which exceed the temperature 
of the cool background by about 3000 K. 

In the magnetic case (Fig. [7]), the shock fronts are virtually 
absent, the divergence of the horizontal velocities typically being 
an order of magnitude smaller than in the non-magnetic simula- 
tion. Apparently, the strong magnetic field suppresses the devel- 
opment of the shock structures. Instead, a small-scale pattern of 
hot filaments appears in the magnetic regions, which is tightly 
connected to vertically orientated vortices and localized up- and 
downflows. The high-temperature features are often associated 
with the shear part of vorticity. 
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Figure 5. Histograms of various quantities in swirling regions for three layers of 210 km width each, centered at heights z = 
-630 km, z = km, and z = 630 km. Solid lines represent the magnetic case, dashed lines the non-magnetic case. The quantities 
considered are swirling period (upper left panel), gas pressure normalized to its horizontal average (upper right), temperature 
normalized to its horizontal average (lower left), and vertical velocity component (lower right). Each histogram is normalized 
by its integral. 



As an example of the flow structures in the magnetic simu- 
lation, Fig. [9 depicts a pair of photospheric vortices in the mag- 
netic case. The velocity streamlines are helically wound with a 
pitch lower than the height of the photosphere. The magnetic 
field is expanding with height and not significantly twisted. This 
is a consequence of the high Alfven speed in the photosphere, 
which reaches about 70 km s" 1 at the upper boundary of the com- 
putational box: any twist of the vertical magnetic flux concen- 
trations rapidly escapes upward in the form of torsional Alfven 
waves. Such processes may drive swirling flows in the chromo- 



sphere ( Wedemeyer-Bohm & Roupp e van der Vo ort 2009 ) and 
in the corona (Zhang & Liu|20lT ). 



3.3. Heating processes 

Figure [K)| shows height profiles of the temporarily and horizon- 
tally averaged temperature and dissipation rates (viscous and 
Ohmic). While the two simulations have very similar tempera- 
tures in the deep and near-surface layers, the average atmosphere 
of the magnetic run shows a much flatter temperature profile 
above z - 200 km. At z = 600 km, the temperature difference be- 
tween the two runs reaches about 600 K. Above this height, the 
non-magnetic atmosphere shows a steep temperature increase, 
most likely due to shock heating. The flat profile in the mag- 
netic case indicates a different heating process that extends over 
a bigger height range and affects also the middle photosphere. 

In order to evaluate the importance of viscous and Ohmic 
dissipation as heating processes, we determined the correspond- 
ing specific dissipation rates, i.e., the change of internal energy 



per unit time and unit mass due to viscous and magnetic dif- 
fusion, respectively. The MURaM simulations analyzed here were 
carried out with a physical magnetic diffusion term in the in- 
duction equation, assuming a constant magnetic diffusivity of 
7] = 1.1 x 10 11 cm 2 s" 1 , the smallest value compatible with the 
given spatial grid resolution. In contrast, the viscosity term in 
the momentum equation effectively acts only near the grid scale. 
It involves artificial diffusivities explicitely depending on the 
local velocity structu re, namely hyperdi ffusivities and shock- 
resolving diffusivities ( Vogle r"et al.|2005| ). For stability reasons, 
there is also magnetic hyperdiffusion near the top and bottom of 
the computational box, but this does not significantly contribute 
to the dissipation rate below z - 700 km. 

The bottom panel of Fig. [TO] shows the height profiles of the 
horizontally averaged specific dissipation rates for both runs. 
The viscous dissipation rate in the non-magnetic case rises 
steeply above z - 600 km, probably related to shock formation 
in this height range. In the magnetic simulation, the dissipation 
is dominated by Ohmic dissipation below z — 500 km, which ex- 
ceeds the viscous dissipation by up to a factor of 5 in the middle 
photosphere. This extended range of heating explains the rela- 
tively flat average temperature profile in the magnetic case. The 
continuous increase of viscous dissipation above z - 200 km is 
partly a result of the assumed height profile for the hyperdiffu- 
sivity, so that it is possibly overestimated. 

Figure [TTJprovides information about the spatial distribution 
of the dissipation rates in the higher layers. In the non-magnetic 
run, the main sources of viscous dissipation are shock fronts. 



5 



temperature [K] 



Moll et al.: Vortices, shocks, and heating in the solar photosphere 



5500 



3 - 



rO.OO 



3 4 
[Mm] 



[Mr, 



swirling strength A ci [s 1 ] 



-■0.07 



-a 



5 - 



/ 



Jo. 00 



-0.35 



3 

[Mn 



3 4 

[Mm] 



vorticity | Vx 



- 2A ri IV 



z = 71 4 km 









0.28 


6 




















0.23 


5 






./ 




0.19 


4 






















0.1 4 


> 3 








r ' 














- 0.09 


2 






















1 0.05 


1 










_ 0.00 








0.05 



3 4 
[Mm] 



3 4 
[Mm] 



Figure 6. Horizontal maps of various quantities at z = 714 km for the non-magnetic case. From left to right, top to bottom: temper- 
ature, vertical velocity, swirling strength, divergence of the horizontal velocity, vorticity, and shear part of vorticity. 



In the magnetic case, the viscous dissipation^] is associated with 
the vortices located in the magnetic flux concentrations while 
the Ohmic dissipation rates are concentrated at the edges of the 
magnetic flux concentrations. 

More detailed and quantitative information on the relation 
between temperature and flow structure in the upper layers is 
obtained from the bivariate (2D) histograms shown in Fig. [12 



2 Since the hyperdiffusivities generally are anisotropic, the associated 
dissipation rates are not necessarily always positive. A few grid cells 
indeed show slightly negative values of the viscous dissipation rate (see 
the small isolated white patches in the middle panel of Fig. [TTJ, but 
these are completely irrelevant for the distribution and the average of 
the dissipation rates. 



The clearest difference between the non-magnetic case (left col- 
umn) and the magnetic case (right column) is in the role of the 
velocity divergence: while there is a strong correlation of nega- 
tive divergence (associated with shocks) with temperature in the 
non-magnetic case, such a relation is almost absent in the mag- 
netic run. In both cases, the temperature correlates with swirling 
strength, but cells with moderate to high swirling strength are 
much more abundant in the magnetic case. In the non-magnetic 
simulation, the correlation with swirling strength is due to the as- 
sociation of swirls with shock fronts. Note that the non-magnetic 
simulation has a much larger proportion of very cool regions 
with temperatures between 2000 K and 3000 K, which are al- 
most absent in the magnetic case. Both simulations show cells 
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Figure 7. Horizontal maps of various quantities at z = 644 km for the magnetic case. From left to right, top to bottom: temperature, 
vertical velocity, swirling strength, divergence of the horizontal velocity, vorticity, shear part of vorticity, strength of the vertical 
magnetic field, and strength of the horizontal magnetic field. 
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with high temperature for all values of the swirling strength, but 
temperatures exceeding 7000 K apparently can only be provided 
by the shock fronts in the non-magnetic case. 
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Figure 8. Horizontal cuts in x-direction through the shock lo- 
cated at (x,y) = (4.5 Mm, 5.7 Mm) in Fig. [6] Upper panel: hor- 
izontal velocity (blue, solid line) and temperature (red, dashed 
line); lower panel: density (blue, solid line) and gas pressure 
(red, dashed line). 



4. Summary and discussion 

The simulation results reveal a striking difference between the 
magnetic and the non-magnetic case in terms of the velocity 
structure in the upper photosphere and around the temperature- 
minimum region. While these layers are dominated by shocks 
in the non-magnetic simulation, the magnetic case is character- 
ized by extended vertical vortices associated with the magnetic 
flux concentrations. Clearly, there must be a transition from the 
shock-dominated non-magnetic regime to the vortex-dominated 
magnetic regime. High-resolution simulations of small-scale dy- 
namo action (Pietarila Grah am et al.||2010| ) indicate that the 
shock pattern is almost unaffected by the low-lying magnetic 
loops and occasionally forming vertical flux concentrations gen- 
erated by the small-scale dynam o. A run simula ting the decay of 
a mixed-polarity field (cf. Cameron et al. ||2 11) shows that the 
suppression of the shock pattern sets in for a mean strength of 
the vertical field component of about 20 G in t he lower photo- 
sphere. This is consistent with the simulations of S chaff enberger 
|et al.| ( [2006] ), who find a clear pattern of shock waves for a mean 
vertical field of 10 G. Consequently, 'quiet' regions with a mean 



unsigned vertical field around 10 G (e.g., Lites et aL|2008| ) are 
probably dominated by shock waves, while the upper photo- 
sphere of more strongly magnetized unipolar regions (in network 
and plage) is dynamically dominated by vortices and shear flows 
associated with the magnetic features. 

In both simulations, substantial local heating takes place: 
shock heating in the top layers of the non-magnetic case and 
a combination of Ohmic heating (in the lower to middle pho- 
tosphere) and viscous heating (in the upper photosphere) in the 
magnetic case. The different height variation of the heating rates 
in the two cases leads to a sharp temperature rise in the upper 
layers of the non-magnetic case, in contrast to a more extended 
temperature enhancement resulting in a flatter temperature gra- 
dient in the magnetic run. 

However, while yielding important qualitative insight in the 
different mechanisms for mechanical heating, we do not con- 
sider the dissipation rates and thus the temperature structure in 
the higher layers provided by the simulations to be quantitatively 
reliable. In the simulation runs considered here, the hyperdiffu- 
sive viscous dissipation rate was arbitrarily enhanced in the up- 
per photosphere for reasons of numerical stability, which affects 
the temperatures. It is also not clear how strongly the dissipa- 
tion rates depend on the spatial grid resolution and on the as- 
sumption of a vertical field at the upper boundary (as opposed 
to matching to a potential field, for instance). Apart from these 
problems, we have to keep in mind that the conditions in the lay- 
ers around the temperature minimum probably deviate from the 
gray LTE radiative transfer used in the simulations. Furthermore, 
a direct comparison between horizontal temperature averages 
(or even averages over surfaces of constant optical depth) from 
simulations with semi-empirical models would be problematic 
owing to the strong horizontal inhomogeneities and the nonlin- 
ear averaging underlying obs erved spectral mean line profiles 
flUitenbroek & Criscuoli|2011| ). 

We intend to continue this study by carrying out new sim- 
ulations with an im proved treatment of the subgrid scales (see 
|Rempel et al.|2009] ) and by performing a detailed investigation 
of the dependence of the dissipation rates on grid resolution and 
on the upper boundary condition for the magnetic field. 
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Moll et al.: Vortices, shocks, and heating in the solar photosphere 






solid: magnetic case 
4000 - dashed: non-magnetic case 




Figure 9. Streamlines (colored) and magnetic field lines (gray) 
of two counter-rotating vortices above the optical surface. The 
vortex pair is located at (x,y) ^ (2.0 Mm, 2.9 Mm) in Fig. [7] 
Color represents the temperature. 



Figure 10. Height profiles of the horizontally averaged temper- 
ature over the full height range (top panel) and above the optical 
surface (middle panel). The bottom panel shows the horizontally 
averaged viscous dissipation rates (solid: magnetic case; dashed: 
non-magnetic case) and the Ohmic dissipation rate in the mag- 
netic case (dotted). 
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Figure 11. Dissipation rates on horizontal cuts through the upper 
photosphere. Top: viscous dissipation in the non-magnetic case; 
middle: viscous dissipation in the magnetic case; bottom: Ohmic 
dissipation. 
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